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Abstract 

Unintegrated parton density functions (uPDFs) of Quantum Chro- 
modynamics (QCD), also known as TMDPDFs, are generally used to 
study details of hadronic final states in high energy lepton-hadron and 
hadron-hadron collisions; while the integrated parton density func- 
tions (PDFs) are used for conventional deep inelastic inclusive pro- 
cesses. The self-similarity based Model of proton structure function 
F2(x,Q 2 ) suggested in recent years are however based on specific 
uPDFs with self-similarity at small x. In this work, we study large x 
limit of such a Model and modify the defining uPDFs to make them 
compatible with theoretical expectations in such limit. Possibility of 
saturation of Froissart bound in this Model is discussed. We also re- 
analyze the PDFs in this approach using its conventional relation with 
TMDPDFs. 
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1 Introduction 



Since nineteen eighties the notion of fractals has found its applicability in high 
energy physics through the self-similar nature of hadron multi-particle pro- 
duction processes [H [2l [3] . Specifically during mid-nineties, James D. Bjorken 
jl] highlighted the fractality of parton cascades leading to the anomalous di- 
mension of the phase space. Similarly in early twenties, fractal characters of 
hadrons had been pursued within a statistical quark model [5] with consid- 
erable success. 

Relevance of these ideas in the contemporary physics of deep inelastic 
scattering (DIS) was first noted by Dremin and Levtchenko [6] in early 
nineties where it was shown that saturation of hadron structure function 
at small x may proceed faster if the highly packed regions (hot spots) of pro- 
ton have fractal structures. However, these ideas received wider attention in 
2002 when Lastovicka [TJ [5] proposed a relevant formalism and a functional 
form of the structure function F 2 (x, Q 2 ) at small x. In recent years, the for- 
malism was further analyzed phenomenologically [9j [TOj [TT] . 

The study of structure functions at small x has become topical in view of 
the high energy colliders like HERA [T2] and LHC []31 EH] where previously 
unexplored small x regime is being reached. In fixed-target DIS experiments 
scaling violations have been observed, i.e. the variation at fixed values of 
Bjorken-x of the structure functions with Q 2 , the squared four-momentum 
transfer between lepton and nucleon. These scaling violations of F2(x,Q 2 ) 
are well described by the DGLAP evolution equations [151 HI E] • The strong 
scaling violations observed at low x are attributed to the high gluon density 
in the proton. In the Quark Parton Model (QPM), F 2 (x,Q 2 ) is the sum 
of the quark and anti-quark x distributions, weighted by the square of the 
electric quark charges, i.e. 

F 2 (x, Q 2 )=xJ2 c?(ft(x, Q 2 ) + Qi{x, Q 2 )) (1) 

i 

The formalism of Ref. [TJ |S] was based on unintegrated parton density 
fi(x, k 2 ), or presently more familiar, Transverse Momentum Dependent Par- 
ton Density Function TMDPDF [18], where k 2 is the parton transverse mo- 
mentum squared. Structure function F 2 (x,Q 2 ) was then obtained by using 
the relation 

Q 2 

q t (x,Q 2 ) = J dklf^k 2 ) (2) 
o 
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As TMDPDF fi(x,kf) is by construction dimensionless, the integrated 
one qi(x, Q 2 ) becomes dimension of inverse area contrary to the common ex- 
pectation. 

One of the main objectives of the present paper is to study if this formal- 
ism can be reformulated by using the standard relation [I§J [2U1 EE] : 




o 



instead of Eq (2) so that the integrated quark density also becomes dimen- 
sionless. 

One of the boundary conditions of the structure function F 2 (x, Q 2 ) is that 
it should vanish as x -> 1 [22j [231 [2j] , i-e. 

limF 2 (a;,Q 2 ) = (4) 

x— >1 

In the present work we will explore if the Fractal Inspired Model of Ref. 
[JJ [8] is compatible with this condition. We will show that this boundary 
condition is not physically possible because the function becomes indetermi- 
nate. We will therefore go back to the original TMDPDF and modify it by 
introducing a new function h(x, k 2 ) in it. 

We also explore the possibility of incorporation of saturation of Froissart- 
Martin limit [25| [26] within this approach as it has attracted attention in 
the recent literature [27J ESI 122] • It is well-known that in the conventional 
QCD evolution equations like DGLAP [T51 Ml E] and BFKL gU\ ED E2J, 
this limit is violated; while in DGLAP approach the small x gluons grow 
faster than any power of lnl/x ~ ln(s/Q 2 ) [22], in BFKL it grows as power 
of l/x [301 ED S21 [33]. 

The study of the Froissart-Martin bound in DIS is more complicated than 
that of the hadron-hadron case because we have the external virtual photon 
mass Q as an additional kinematic variable. A remarkable modification of 
the bound in DIS analysis has been suggested by Gribov [31]. Later Gots- 
man, Levin and Maor |35j proposed a generalization of the Gribov formula. 
We will comment on the implication of such generalization in the TMDPDFs. 

In the Ref. [JJ [8] the basic TMDPDFs have only one hard scale k 2 (as 
in BFKL approach) but in general, it may have two hard scales: k 2 and the 
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photon virtuality Q 2 [20] as in the CCFM approach EZ] • 



In this paper, we will also generalize the formalism of Ref. [7, 8] with 
TMDPDFs having two hard scales. 

2 Formalism 

2.1 Self-similarity based TMDPDFs with One Hard 
Scale 

One of the most intriguing problems in QCD is the growth of the cross- 
sections for hadronic interactions with energy. Conventionally deep inelastic 
lepton-hadron scattering is described in terms of scale-dependent parton dis- 
tributions xg(x,k 2 ) and xq(x,k 2 ). These distributions correspond to the 
density of partons in the proton with longitudinal momentum fraction x and 
parton transverse momentum k 2 [T?] . 

The self-similarity based Model of the nucleon structure function pro- 
posed in Ref. [7J [8] has been designed to be valid at small Bjorken-x. The 
formalism described in these references is based on the imposition of self- 
similarity constraints to the dimensionless TMD quark density fi(x, kf) and 
relate it to the integrated density. In other words, using magnification factors 
1 , / k 2 



x V Q 



— and ( 1 + -—^ J , an unintegrated quark density (TMD) is given as: 

Dl log \ log i 1 + Q 2 ) + ° 2 log I + Ds log { 1+ Ql) + Dl ° 

(5) 

where i denotes a quark flavor. Here, D2 and D3 are the fractal parame- 
ters; Di is the dimensional correlation relating the two magnification factors; 
while .Dq is the normalization constant. Conventional integrated quark den- 
sities [PDF] qi(x,Q 2 ) are defined as sum over all contributions with quark 
virtualities smaller than that of the photon probe Q 2 . Thus fi(x, k 2 ) has to 
be integrated over k 2 (Eq (2)) to obtain qi(x,Q 2 ). 

The following analytical parameterization of a quark density is obtained 
by using Eq (2). 

. 2 e^Qlx-"" ( -Wi+g) /, , Q 2 \ D ' +1 ,, ., . 
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Now using Eq (6) in Eq (1), the expression for proton structure function is 
obtained as: 



where D = Sj-Dq . 

As we have seen above, Lastovicka's Model starts with self-similar TMD- 
PDFs. Imposing large x limit, i.e. at x — >■ 1, Eq (7) becomes: 

F 2 (1,Q 2 ) = ^(( i-:-^) -1) 



Let 

e D oQ 2 1 + D 3 




If F 2 (1,Q 2 ) = 0, then we get _D 3 = —1 and subsequently F 2 (x,Q 2 ) becomes 
indeterminate (using Eq (9)). It implies that the Model does not have correct 
x — > 1 behavior. This calls for generalization of PDFs to add a new term 
in TMDPDFs. A plausible way to incorporate correct x — > 1 behavior for 
TMDPDF is to introduce an additional term h(x, kf) which does not have 
fractal basis. 

Now, starting with the TMDPDF, we have: 



logMx,k?) = D 1 log 1 log (l + ^+D 2 log ± + £> 3 log + ||J + 
Dl + ]ogh(x,k*) (10) 



i.e. 



/.<*,*>=(;) ^g) + e D °h(x,k1) (11) 



Thus the integrated quark density is given as: 



ft(x,Q 2 ) = y d*?^) (l + ^J e D oh(x,kf) (12) 
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In the simplest case, consider h(x, k 2 ) to be only ^-dependent, i.e. h(x, kf) w 
h(x). 

Then, the improved model is 

F*(x,Q 2 ) = F 2 {x,Q 2 ).h{x); where 0<x<l (13) 

where F 2 (x, Q 2 ) is the original Lastovicka's model (Eq (7)). 

Neglecting kf dependence of h(x, k 2 ), the most economical way of having 
correct large x behavior is h(x,kf) pa h(x) ~ (1 — x) a ,for a > 0, but in 
general the integration is to be carried out with its explicit kf dependence in 
Eq (12). Thus large x effect should be introduced as an additional component 
of TMDPDF to make the PDF compatible with the theoretical expectation, 
i.e. F 2 (x, Q 2 ) -»■ as x -> 1 (221 1231 I2H- 



2.2 Self-similarity based TMDPDFs with Two Hard 
Scales 

In general, TMDPDF depends on two hard scales, parton transverse mo- 
mentum kf and the scale Q 2 of the probe [201 EH GS] • In Fractal Inspired 
Model (Lastovicka's), it will therefore have three magnification factors, viz. 
I { k 2 \ { Q 2 \ 

~i ( 1 + "tt J an d I 1 + 7^7 I ; and four correlation terms. The generalization 

x V K) \ QoJ 

of TMDs of this Model is therefore given as: 

2 l ^ l 2\ 



logfi(x,k t ,Q 



1 , ( . kr\ „ . 1 .... ( . kf 



Alog-logl 1 + ^|J + J D 2 log- + J D 3 log(l + - 7 i ) 



D * log{1 + W +D5los x log { 1 + W + 

D s \ogh(x,kf,Q 2 )+Dl (14) 



where h(x, kf, Q 2 )accounts for large x behavior. 



In the above Eq (14), there are altogether eight parameters, viz.D 2 , D 3 , 
(Fractal dimensions) ;Di, D 5 , D e , D 7 (Correlation parameters) and D l Q 
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{Normalization constant); as well as one unknown function h(x,k 2 ,Q 2 ). 

While D\, -D5, Dq, Dj are the correlation parameters correlating the 
different magnification factors, h(x, k 2 , Q 2 ) is the additional function that 
cannot be constrained by the notion of fractal geometry. 
Now following the same steps as above, Eq (14) can be rewritten in the 
following way: 

fi(x, k 2 ,Q 2 ) = ^.^,^Q 2 )(i) D2 (l + |) D3 (l + g) D4 . 

i\ Dilog ( 1+ !) m D5log ( 1+ €> ( . k 2 \ D6 H 1+ ^) 



X J \x 

1^ D710S ( 1+ I) 10g ( 1+ C) 

The integrated quark density is then given as: 

Q 



where 



1 + 



(15) 



^^^log^l+^j ^ Q2 ^D 4 



<h(x,Q 2 ) = e u '°[-) + (1G) 



J(a;,Q 2 ) = | dk 2 h(x,klQ 2 ) 



i\^ io K 1+ i) +D7los ( 1+ l) Iog ( 1+ ^) 



, 2 x D 3 +D 6 log l+« 

1 + " < 17 > 

Eq (16) is the most general form of PDF which contains two cut-off param- 
eters: /c 2 and Ql . 

In a special case, imposing the limit _D 4 , D 5 , D 6 , D 7 = and -D 8 = (i.e. 
h{x, k 2 , Q 2 ) = 1 [see Eq (14)]), the above integral of Eq (17) becomes: 



l\ Dll °s( 1+ f 



I(x,Q 2 ) = j dk 2 ^-j * U/ { 1+ '^J ( 18 ) 


which is the original Lastovicka's Integral to obtain Eq (6). Now let us 
proceed to evaluate the integral of Eq (17), in more general case, with the 
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following observation: 

The exponent of — in Eq (17) can be re-expressed as 

x 

log (l + |) D x with D x = { D 1 + D 7 log (l + gjj) } . 

Similarly, the exponent of ^1 + ^ in Eq (17) is D 3 = D 3 +D 6 log (l + §r) ■ 

Without explicit form of h(x,k 2 ,Q 2 ), Eq (17) cannot be evaluated ana- 
lytically. 

In the limit h(x, k 2 , Q 2 ) = 1, Eq (17) then becomes: 

Q 2 





which is again similar to Eq (18), i.e. Lastovicka's Integral, except that in 
this case D x is replaced by D\ and D 3 by D 3 . Eq (19) implies 



Now substituting this integral from Eq (20) in Eq (16), we obtain: 



qi(x,Q 2 ) = 




The expression for the proton structure function Fl(x,Q 2 ) can thus be ob- 
tained using Eq (21) in Eq (1): 



Fi(x,Q 2 ) = 
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where D\ and D 3 are independent of k 2 . 

As in Eq (7), imposing boundary condition (Eq (4)) at x — > 1, Eq (22) 
yields D 3 = — 1 and subsequently F 2 7 (x,Q 2 ) becomes indeterminate. 

Thus the correct x — > 1 limit of the structure function (Eq (22)) cannot 
be incorporated as in the original one. Therefore the most economic way to 
incorporate large x behavior is by putting h(x, k 2 ) ^ 1 but of the factorizable 
form h(k 2 )(l—x) a so that it vanishes as x — > 1. As in the original Model it can 
be seen that x — > 1 limit of Eq (22) with h(x, k 2 ) ~ h(k 2 )(l — x) a , for a > 0, 
exists. As noted earlier (with Eq (13)), we have neglected explicit k 2 (and 
Q 2 as well) dependence of h(x, k 2 , Q 2 ) in the most economical version of the 
model. However, as noted above, the integration of Eq (17) can be carried 
out only with its explicit k 2 dependence. 

2.3 Problem of Dimensionality 

Exact relation between TMDPDF and PDF is still an active area of study {3"8l 
[391110]. Formally, the integration over the transverse component of the parton 
momentum is expected to yield the PDF (Eq (2)). However, the self-silmilar 
TMDPDF (Eq (5), Eq (10) and Eq (14)) are by construction dimensionless in 
the formalism. As a result, the resulting PDF as well as the proton structure 
function possess dimension of inverse area. Literally speaking, PDF has 
become dimensional while TMDPDF is dimensionless. Phenomenologically, 
it matters not much but in principle, it is questionable. If instead we consider 
Eq (3), k 2 in the denominator plays the role of dimensionality check. But even 
if dimensionality problem is solved, there arises a singularity at k 2 = 0. This 
can be overcome by introducing a minimum transverse cut-off momentum. 
The crucial point is that in order to make the PDF finite and dimensionless 
in the formalism, we need to have a cut-off over k\ , say k 2 . 
Thus, the dimensionless Lastovicka's integral can be written as: 

=/itG) (23) 

So the dimensionless PDF in the original Lastovicka's Model becomes: 

q :(x,Q 2 ) = e D ^(-Y 2 1(x,Q 2 ) (24) 



9 



and the corresponding structure function becomes: 

F;(x,Q 2 ) = e D °(±y 2 1 1(x,Q 2 ) (25) 

to be compared with Eq (6) and Eq (7) respectively. 

The formalism can be generalized to TMDPDF with two hard scales 
as well. In this case, dimensionless PDF in the two-scale generalized Las- 
tovicka's Model is given as: 

/-, \ D 2 +D 5 log(l+^ ) , n ^\ Di 

where 



and consequently, the structure function has the form: 



/ , \ D 2 +D 5 log 1+ « -1 / n2 v D 4 

F 2 (x,Q 2 ) = e^Qj V Q ° ; + I(x,Q 

to be compared with Eq (21) and Eq (22) respectively. 



(2f 



2.4 Froissart Bound and Self- Similarity 

Froissart bound [251 ES] implies that the total cross-section (or equivalently 



structure functions) cannot rise faster than the logarithmic growth of In 2 




Such a slow logarithmic growth is however not compatible with the notion 
of fractal dimension as occurred in the original Fractal Inspired Model [7J |8]. 
The reason is that self-similar models allow only power law growth and not 
the logarithmic one. This statement can be explained in a more explicit way 
in the following manner. 



(i) Power law growth 



Let N(R) be the number of self-similar objects, R{N) be the magnifica- 
tion factor and D the exponent of the power law behavior, then the relation 



10 



is given as: 

N(R) = cR(N)- D (29) 
where c is the proportionality constant. Eq (29) leads to 

fls3 lim W> (30) 

N^oo,R-+0hg(l/R(N)) 

Eq (30) shows that D is ultimately independent of N(R) and R{N) in the 
limit N — > oo and R — > 0. It implies that 17 is universal property of a self- 
similar object to be identified as its fractal dimension [%!] . 

(ii) Logarithmic growth 

If, on the other hand, number of self-similar objects is assumed to grow 
logarithmically, i.e. 

NW = <w (31) 

one obtains 

In N(R) , . 

53 -^m < 32 > 

Eq (32) implies that limiting D is not N(R) or i?(iV)-independent and hence 
it is not an intrinsic property of the object. So it cannot be identified as 
fractal dimension. 

Thus the logarithmic growth and hence the Froissart behavior cannot be 
described in a self-similar model. 



2.5 Froissart Bound and TMDPDFs 

Understanding the behavior of high energy hadron reactions from a funda- 
mental perspective within QCD is an important goal of high energy physics. 
It was already observed that the total hadronic cross-section grows with the 
centre of mass energy (s/s). A QCD based explanation for this rising be- 
havior was proposed by Gaisser and Halzen [12]: the cross-section would 
grow because partons would begin playing a dominant role in the hadronic 
reactions. At high energy, the parton distribution, especially the gluon dis- 
tribution, grows very rapidly and thus leads to rising cross-section. It is 
expected that at very high energies the hadronic cross-sections satisfy the 
Froissart bound, which is a well-established property of the strong interac- 
tions [2S]. It has been shown that the Froissart bound is saturated at very 
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high energies [4"3"| HI] in jp, n^p and pp & pp scattering. Saturation of the 
Froissart bound refers to an energy dependence of the total cross-section ris- 
ing no more rapidly than In s [27 \. 

It is already known that the BFKL evolution equation [301 13T1 132] sums 
leading logarithms ln(l/x). The gluon distribution function, governed by this 
equation, is found to increase at small x and it has been confirmed by the re- 
cent HERA data of the proton structure function F 2 (x, Q 2 ) which is involved 
in DIS [15], 146] . The increase is however power-like, such that F 2 (x,Q 2 ) and 
the DIS cross-section rise as power of (^), i.e., as power of s for x = Q 2 /s. 
This behavior does not satisfy the Froissart bound and thus violates unitarity 
[47] . The BFKL equation cannot thus be the final theory for small x physics. 

This bound [251 ES] derived from analyticity and unitarity demands that 
F 2 (x, Q 2 ) grow no more rapidly than In 2 f — ] at very small x for all values 

w 

of Q 2 . Therefore, for fixed Q 2 and small x, if the condition 

F 2 (x,Q 2 )h(x,k 2 ) < ln2 (^) (33) 

is satisfied, then the Model satisfies the Froissart bound with F 2 (x,Q 2 ) and 
h(x, k 2 ) as defined in Eq (7) and Eq (11). 

Consider F 2 (x, Q 2 ) at Q 2 = Qq and ultra small x. Then, -Dilog >> 
D 3 , 1 and therefore Eq (7) yields 

/l\D 2 +D! In2-1 

F 2 (x,Q 2 )* ln(1) (34) 

Using Eq (34) in Eq (33), we have: 

In 3 (I) 

h( X ,Qo)~ / 1 xD 2 +Diln2-l ( 35 ) 
\x) 

for a given k 2 . 

Eq (35) gives the specific limiting behavior of h(x, Q 2 ) for Froissart bound 
compatibility of Eq (7) at Q 2 = Qq. 
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Froissart compatibility condition (Eq (35)) can be generalized to any Q 2 
for ultra small x. 

In 3 (±) 

h(x, Q 2 ) < j^- (36) 

In Ref. [35], Gotsman, Levin and Maor had improved Froissart-Martin bound 
for ultra small x at fixed Q 2 for virtual photon-proton cross- sect ion: 



Fii^Q 2 ) < f ln B/2 1 

1 — X ~ \ % 



*(r\p)*^P^< In 5 / 2 - (37) 



In that case the corresponding improved Froissart compatibility condition on 
the TMDPDF factor h(x, k 2 ) will be 



In 3 (I) 



for any given k 2 at Q 2 = Q\- 

2.6 Graphical Representation of Self-similar TMDPDFs 

Specific models of TMDPDFs available in recent literature are summarized 
in Ref. [48J. Here we analyze the ones of Ref. [7J |8] alone. The TMDPDFs 
fi(x,k 2 ) of the original Lastovicka's model (Eq (5)) contain the parameters 
Dq, Di, D 2 , D 3 and Qq. These are determined from HERA data in the range: 

0.045 < Q 2 < 150 GeV 2 
and 6.2 x 10~ 7 < x < 0.2 (39) 

where 



D = 


0.339 ±0.145 




D 1 = 


0.073 ±0.001 




D 2 = 


1.013 ±0.01 




D 3 = 


-1.287 ±0.01 




Ql = 


0.062 ± 0.01 GeV 2 


(40) 


g n f = 


4, Di w 0.085 (taking 


,D = Z i Di ) =n f Di ) ), 



fi(x, k 2 ) is given by Eq (5), viz. 



r r > » i i\ D2+Dll i 1+ i)( 1+ AT 3 

f i (x,k 2 t ) = e D ° - (41) 



x 
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In Figure 1, fi(x,kf) vs kf is shown for representative values of (i) x = 
10" 4 (ii) x = 10- 3 (iii) x = 1(T 2 and (iv) x = 0.1. 
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Figure 2: fi(x, k%) vs x 

3 Conclusion 

In this paper we have generalized the method of parameterization of the 
self-similar TMDs of Ref. [7J [8] to two hard scales. We have shown that 
the redefinition of the TMDs is necessary to make them compatible with the 
theoretical expectation in large x limit. We have then studied the compati- 
bility of saturation of the Froissart-Martin bound within this approach and 
obtained necessary constraints on TMDPDFs for it. Relation between the 
TMDPDFs and PDFs within such Fractal Inspired Model is also re-examined 
in the context of their relative dimensions. 

Let us now compare the present approach with other models based on 
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self-similarity and fractals. In the statistical quark model of Ref. [5], the 
proton possesses a fractal structure of fractal dimension 9/2, topological di- 
mension 3 and anomalous dimension 3/2. The present approach falls short 
of such specific predictions as it has inbuilt multifractal characters. As noted 
earlier, in the original version of the Fractal Inspired Model [TJ |8], one has 
two fractal dimensions D 2 , D 3 and one correlation parameter D\ fitted from 
HERA data. In Ref. [49] . the monofractal limit of the model [7J |S] was 
indeed studied. Its analysis however indicated that only in a limited x — Q 2 
range (Q 2 < 12 GeV 2 ), such notion of monofractality makes sense. In that 
range, dimensional correlation (Di) vanishes and the proton possesses frac- 
tality (D « 1.107) close to Koch curve (D « 1.26). 

In Ref. [B] , the proton is viewed as a collection of smaller colorless sources, 
filling its interior. The density of such sources has fractal structures which 
gives rise to faster saturation of the proton structure function at small x. 

It is also instructive to recall the physical interpretation of fractal di- 
mension of proton since the notion is rather recent in its description. As is 
well-known [50J, the fractal dimension measures the way in which distribu- 
tion of points fill a geometric space on the average. Extending the notion to 
TMDPDF, fractal dimension tells how densely small x partons fill the proton 
in self-similar way on the average. It is also interesting to note that in the 
model of Ref. |T, 8], TMDPDF takes the simplest form 



in the limit D = D 2 >> D\, -D3.R1 such fractal dimension is essen- 

tially close to x-slope [51] or Pomeron intercept [521 [531 154"] - 

To conclude the problem of self-similarity in the gluon cascades is an 
interesting one and worth investigating. In this paper, we have used fractality 
only as a tool to provide parameterization for the TMDPDFs having one and 
two hard scales. These result in a set of parameters to be fitted from data, 
as was the case with the original model of Ref. [7J |8]. However, in spite 
of this inherent limitation, it has been shown how large x effects and high 
density QCD phenomenon like Froissart saturation should appear in such 
fractal inspired TMDPDFs. Naturally the most interesting study would be 
to calculate the fractal parameters from the first principles, say from a self- 
similar Lagrangian [55]. Interesting though, the solution appears to be a 
formidable task, beyond the reach of the present authors. 




(42) 
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